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Bell’s theorem sets a boundary between the 
classical and quantum realms^, by providing a 
strict proof of the existence of entangled quan¬ 
tum states with no classical counterpart. An ex¬ 
perimental violation of Bell’s inequality demands 
simultaneously high fidelities in the preparation, 
manipulation and measurement of multipartite 
quantum entangled states. For this reason the 
Bell signal has been tagged as a single-number 
benchmark for the performance of quantum com¬ 
puting devices^“^. Here we demonstrate deter¬ 
ministic, on-demand generation of two-qubit en¬ 
tangled states of the electron and the nuclear spin 
of a single phosphorus atom embedded in a sili¬ 
con nanoelectronic device^. By sequentially read¬ 
ing the electron and the nucleus, we show that 
these entangled states violate the Bell/CHSH^ in¬ 
equality with a Bell signal of 2.50(10). An even 
higher value of 2.70(9) is obtained by mapping 
the parity of the two-qubit state onto the nuclear 
spin, which allows for high-fidelity quantum non¬ 
demolition measurement (QND)^ of the parity. 
Furthermore, we complement the Bell inequality 
entanglement witness with full two-qubit state to¬ 
mography exploiting QND measurement, which 
reveals that our prepared states match the tar¬ 
get maximally entangled Bell states with >96% fi¬ 
delity. These experiments demonstrate complete 
control of the two-qubit Hilbert space of a phos¬ 
phorus atom, and show that this system is able 
to maintain its simultaneously high initialization, 
manipulation and measurement fidelities past the 
single-qubit regime. 

Bell’s theorem provides a boundary to the strength 
of correlation that a pair of quantum two-level systems 
(qubits) can display, under the assumption that physical 
systems cannot be instantly affected by distant objects 
(“locality”) and that their properties exist before they 
are observed (“realism”). Certain quantum entangled 
states are predicted to violate Bell’s theorem, and there¬ 
fore invalidate local realistic interpretations of quantum 
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mechanics^. The most profound implications of Bell’s 
theorem arise when observing entangled pairs of parti¬ 
cles that are separated in space-time, such as photons 
travelling at the speed of light in different directions^’ 
Even for physical systems where space-time separation is 
not achievable, the Bell test carries an important mes¬ 
sage when placed in the context of using such systems 
for quantum information processings^. It is in princi¬ 
ple possible to produce and observe quantum entangle¬ 
ment without violating Bell’s inequality, as is the case 
in a wide range of experiments, particularly in the solid 
stateS^~S7 However, a Bell’s inequality violation consti¬ 
tutes an even stronger measure of the ability to faithfully 
produce, control and read out non-classical states of pairs 
of qubits^”^. This maps directly onto the ability to per¬ 
form high-fidelity entangling operations'^, which, along 
with single-qubit operations can fully access the full two- 
qubit Hilbert space. 

Experimental access to Bell’s theorem takes the form 
of the Clauser-Horne-Shimony-Holt (CHSH) inequality^. 
It involves the joint measurement of a two-qubit system 
along measurement axes a and f3. The binary measure¬ 
ment outcomes (0/1) produce a correlation 

E{a, p) = Poo + Pii — Poi — Pio (1) 

where P is the probability of detecting the subscripted 
measurement outcome. In its standard form, the inequal¬ 
ity is tested by measuring each qubit along two axes, a, a' 
and P, P' respectively, and extracting E in all four possi¬ 
ble combination of axes. If the measurement setup does 
not allow for physically rotating the measurement axes, 
equivalent outcomes are obtained by rotating the qubit 
prior to a measurement along a fixed axis. The Bell signal 
is then 

5 = E{a, p) + E{a\p) + Eia, p') - E{a\ p'). (2) 

Bell’s theorem states that, within local realistic theo¬ 
ries, \S\ < 2. Conversely, quantum mechanics predicts 
Smax = 2v^ for a maximally entangled state and an ap¬ 
propriately chosen sets of axes. 

In the present experiment we use as qubits the elec¬ 
tron (I / t) and the nuclear (J/ / spins of a single 
substitutional donor, implanted^^ in an isotopically 
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Figure 1 | Device operation and state preparation protocols, a. Coloured scanning electron microscope 
image of a device similar to the one used in the experiments. The gates coloured in red make up the single electron 
transistor (SET), used to perform sensitive charge sensing for electron spin readout; the donor gates (DG) coloured 
in blue, control the donor potential; the gate coloured in yellow is a broadband antenna used to drive the nuclear 
magnetic resonance and electron spin resonance (NMR/ESR) transitions. Inset shows a schematic four-level energy 
diagram of the two-spin electron-nuclear system under a static magnetic field. The colour coding of the three 
accessible transitions (mwI, , rf) will be used in the pulse sequences throughout the text and figures, 
b. (Bottom) Electron ||) initialization is achieved by setting the donor potential such that > //set > Since 
jeBo ^ kBT, where /cb is the Boltzmann constant and T « 100 mK is the temperature of the electron reservoir, the 
ionized donor is predominantly neutralized by a |^) electron. (Top) After initialization, is increased to lower the 
donor potential to <C jJiSET^ preventing the electron from escaping the donor during the control phase, c. The 

nuclear spin is initialized by applying the sequence tt If the nuclear spin is jf^), tt flips the electron 

spin to It), after which ttrf is off-resonance and leaves the nuclear spin remains in the target state jf^). If the nuclear 
spin is |J|), TT is off-resonance, allowing the subsequent ttrf to flip the nuclear spin to \i[). A final electron 
\downarrow) initialization phase completes the |tlT) preparation. After the two-spin initialization sequence is 
performed, one of the maximally entangled states (DQC or ZQC) is prepared, d. Electron Rabi oscillations 
highlighting the improvement in visibility with the implementation of selective initialization (see main text) and 
QND measurement protocols. 


enriched epilayer^’^^. The qubits are coupled by 
the hyperflne interaction A ^ 96.9 MHz (shifted from 
the bulk value of 117 MHz due to the strong electric 
fields in the nanostructure^^). A static magnetic field 
Bq ^ 1.55 T induces a Zeeman splitting on the electron 
(7eHo with 7e ~ 27.97 GHz/T) and nuclear (7 nHo with 
7 n ~ 17.23 MHz/T) spins, resulting in the two-qubit en¬ 
ergy levels diagram depicted in the inset of Eigure la. 
Manipulation of the full two-qubit Hilbert space requires 


access to at least three eigenstate transitions. Our ex¬ 
perimental setup includes three microwave sources which 
deliver signals to a nanoscale broadband antenna^^, pro¬ 
viding access to both electron spin resonance (ESR) tran¬ 
sitions (z^Mwi, = 7e^o T ^/2) and one nuclear magnetic 
resonance (NMR) transition (z^rf = 7n^o + ^/2). Goher- 
ent qubit rotations are denoted as ag where a is the rota¬ 
tion angle in radians and the subscript identifies the tran¬ 
sition. The donor is located in the vicinity (~ 25 nm) 
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of the large electron island of a single-electron transistor 
(SET), formed by biasing a metal-oxide-semiconductor 
(MOS) gate stack above a Si02 layer. An additional 
set of gates (DG) controls the (spin-dependent) electro¬ 
chemical potential of the donor, with respect to the 
SET potential /isET (Eigure la)^. The SET island acts as 
a quasi-continuum electron reservoir for spin-dependent 
electron tunnelling to and from the donor, which is the 
key step for single-shot readout and electrical initializa¬ 
tion of the electron qubit state^^. 

On-demand entanglement requires the two-spin system 
to be initialized to a known state. Eigure Ib-c explains 
in detail the initialization procedure. The fidelity is lim¬ 
ited by electron initialization errors. The probability of 
erroneous |t) electron initialization is proportional to the 
thermal population of electrons in the SET island at the 
potential Similarly, there is a small probability for a 
1^) electron to tunnel back to the SET, inversely propor¬ 
tional to the density of states in the SET island at 
Additionally, the finite ratio between the electron tunnel 
rate and the initialization time can be another source of 
error. If the donor is in the ionised state at the time the 
potential is lowered, the loading of the electron occurs 
while both spin states are equally accessible, resulting 
in a random initial electron spin state. To improve the 
initialization fidelity we monitor the SET current dur¬ 
ing the \D initialization phase, and discard the traces 
where the donor is found ionised at the end of that phase 
(see Extended Data Eigure 1). This selective initializa¬ 
tion protocol improves the two-spin initialization fidelity 
from 94.3(7)% to - 96.5(7)%. 

We then apply the two-pulse sequence 7r/2RF:7rMwi or 
7r/2j^p:7r to the state to prepare a maximally en¬ 
tangled Bell state. By selecting MWl or for the sec¬ 
ond pulse, we can prepare the zero quantum coherence 
(ZQC) |T) = ^(|^J|) + lift)) or the double quantum co¬ 
herence (DQC) |4>) = ;^(|tlT) + |iJ|)) states, respectively. 
The entire initialization and entangled state preparation 
sequence is depicted in Eigure Ic. 

Characterization of the prepared entangled state re¬ 
quires high-fidelity single-shot projective measurements. 
The electron spin readout, based upon spin-dependent 
tunnelling from the donor to the SET island^^, is de¬ 
structive (the electron is lost after readout). Its fidelity, 
here up to 97%^, is limited by the thermal broadening of 
the reservoir. In contrast, the nuclear qubit is measured 
through a quantum nondemolition (QND) method^ by 
mapping its spin state onto that of the electron, through 
an electron-nuclear CNOT operation (e.g. by applying 
TT ). The QND readout allows for repetitive measure¬ 
ment of the nuclear spin state, which reaches a fidelity 
> 99.9% with 30 repetitions^. With access to the full 
two-qubit Hilbert space, any observable in the four-level 
system can be mapped to the nuclear spin in order to 
take advantage of the high-fidelity QND measurement. 


To illustrate this point, we perform coherent Rabi ro¬ 
tations of the electron spin qubit. Their visibility is a 
nontrivial function of both measurement and initializa¬ 
tion fidelities'^. Eirst we perform the measurement by 
direct electron readout, then by mapping the electron 
state onto the nucleus and performing QND readout on 
it. By further applying the initialization protocol de¬ 
scribed above, the overall Rabi visibility improves from 
91% to 96% (Eigure Id), which means we have removed 
over half of the state preparation and measurement errors 
in the system. 

The two-qubit system is now optimized to attempt 
CHSH inequality violations. The single qubit rotations 
required by the CHSH protocol are achieved through 
sequential rotations at each of the qubit’s transitions. 
To perform a (3 rotation on the electron independent of 
the nuclear state, we use the sequence /3 mwi-/^ , since 

both electron frequencies (mwI and ) are simulta¬ 
neously addressable. Conversely, for the nuclear spin we 
can only apply RF within the ||) subspace. Therefore 
a a rotation on the nuclear qubit is performed through 
Q^RF-TTMwi-Tr 'cyrf^ where the ttmwi-tt pulses bring the 
It) manifold into the RE-addressable |t) subspace. 

A maximum Bell signal is achieved when the angle be¬ 
tween the two projection axes on the same qubit (e.g. 
a and a') is 7r/2 and the relative shift between the two 
sets of axes is 7r/4. The correlation E is obtained by in¬ 
dependently measuring the individual qubits, through a 
single-shot readout of the electron spin immediately fol¬ 
lowed by a nuclear spin measurement. By repeating the 
sequence 300 times and logging each set of electron and 
nuclear spin results, we can construct the set of probabili¬ 
ties for each eigenstate {7%^, compute 

E from Equation 1. 

Alternatively, we can obtain P in a single measurement 
by first mapping the parity of the Bell states (even parity 
for DQC, odd parity for ZQC) to the nuclear spin state, 
and then performing a QND measurement on it (See Ex¬ 
tended Data Eigure 2a). This is done by applying a ttrf 
pulse before a nuclear measurement, resulting in with 
probability P-^^ + P^^ and |J|) with probability 
Therefore, the correlation (Eq. 1) becomes P = 1 — 2P^. 
The complete sequence for the CHSH experiment—along 
with our chosen set of projected axes—is shown in Eig¬ 
ure 2a. 

The histograms in Eigure 2b show clear violations of 
Bell’s inequality for both DQC and ZQC states, using 
both standard and QND measurement protocols (see 
numbers in Table 1). We obtain a maximum mean Bell 
signal of 2.70(9) for the ZQC state using QND readout, 
with the error taken as the standard deviation of the 
histogram data. Each histogram is constructed from a 
compilation of ^ 120, 000 single-shot measurements (see 
Eigure 2 and Extended Data Eigure 2b for the break¬ 
down), taken over a period of ^ 11 hours, highlighting 
the long-term stability of this system. The on-demand 


4 



Figure 2 | Bell’s inequality violation, a. Protocol for performing the CHSH experiment in the 
electron-nucleus system. After preparation of an electron-nuclear entangled state, each qubit is rotated to obtain the 
desired combination of projection axes. Then, the correlation is extracted by measuring the qubits either directly or 
via mapping the parity onto the nuclear qubit, and performing a single QND readout on it. A correlation value 
E{a^ /3) is calculated after 300 repetitions of this sequence, b. Histograms of Bell signal values obtained for the ZQC 
(left) and DQC (right) states. After obtaining E for each of the 4 combinations of qubit rotations, a bell signal S is 
calculated (Equation 2). The histograms in the figure are constructed from ^ 100 Bell signal values each. We 
assume a normal distribution to calculate the mean and standard deviation of S displayed in Table 1. 


entangled state preparation and single-shot readout tech¬ 
niques adopted here also avoid the detection loophole^^. 

In addition to the Bell test, we have experimentally 
mapped out the density matrix of the entangled states 
by adopting a tomography method first implemented by 
Mehring et al.^^, where the off-diagonal matrix elements 
of an electron-nuclear state are extracted from measure¬ 
ments of a single observable. We apply geometric phase 
operations to the input state and increment the geomet¬ 
ric phase of the electron and nuclear spins relative to 
each other, through the sequence , where (j) and 

0 are the pulse rotation axes relative to the initial state 
preparation pulse. We sweep the phases by increments 
A0, A6> and map the desired off-diagonal coherence to 
the nuclear observable. The mapping is a pulse sequence 
that depends on the particular coherence being measured 
(see Extended Data Eigure 3). The measured as a 
function of phase increment reveals oscillations with am¬ 
plitude proportional to the off-diagonal density matrix 
element (Eigure 3a). The frequency of the oscillations 
is a linear function of A(/) and A^, which can be cho¬ 
sen such that the recovered oscillations for each coher¬ 
ence have distinct frequencies (Eigure 3b). The diago¬ 
nal elements of the density matrix are obtained from the 


offsets of the tomography signals. Extended Data Eig¬ 
ure 3 provides details on the pulse sequences, frequencies 
and offsets for each of the coherences of the system. The 
main advantages of this tomography method are twofold: 
first, the extracted coherence amplitudes are strict lower 
bounds and hence provide a conservative entanglement 
estimate; second, it only requires measurement of one of 
the qubits, which allows us to exploit the much higher 
fidelity of the nuclear qubit QND measurement. Addi¬ 
tionally, this technique reduces false positive signal con¬ 
tributions which could arise from pulse errors. Each co¬ 
herence evolves according to a distinct frequency, and 
when measuring a particular coherence only its particu¬ 
lar frequency numerically contributes to the tomographic 
result; all other frequency contributions are due to pulse 
errors and are discarded (see Extended Data Eigure 4). 

The density matrices (Eigure 3c, see Extended Data 
Eigure 4 for numerical matrix) extracted using this pro¬ 
tocol show remarkable fidelities E = Pi ^/EP 2 -\/E^ > 
96% when compared to the ideal states. Two addi¬ 
tional density matrix tests, the positive partial transpose 
(PPT)^^’^^ and the concurrence (C)^^, both confirm our 
measured states to be highly entangled. Table 1 shows all 
of the entanglement benchmarks obtained from both the 
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Figure 3 | Density matrix tomography, a. Example of state tomography through application of geometric 
phases. After preparing the ZQC state, we apply the geometric phase map, and then move the |tJ|) (ilT| coherence 
onto the nuclear spin observable, which is then measured in QND mode. The oscillations in the |J|) probability 
arise from the increasing phase accumulated by the |tJ|) coherence through the tomography protocol, and the 
visibility of the oscillations provides a lower bound on the amplitude of the |tJ|) density matrix element, 
b. Fourier transforms of tomography protocol signals obtained after preparing and mapping each of the off-diagonal 
coherences (protocols shown in Extended figure 3). With carefully chosen phase increments (A0 = — and 
AO = — ^tt) we can obtain well separated peaks for each coherence, c. Density matrices extracted for the DQC 
(Left) and ZQC (Right) states, expressed in polar coordinates where an off-diagonal element’s phase is indicated by 
its colour. We have added a global phase correction to the matrices such that the maximum off-diagonal elements 
have zero phase. 


Bell test measurements and density matrix tomography. 


The deterministic preparation and detection of quan¬ 
tum entangled states that violate Bell’s inequality with 
S = 2.70(9) provides a striking demonstration of the po¬ 
tential of donor-based qubits in silicon for quantum in¬ 
formation processing. In particular, it highlights the role 
of the nuclear spin in expanding the size of the available 
Hilbert space and in providing a pathway to high-fidelity 
QND measurements for all observables. Future exper¬ 
iments will further explore the role of nuclear spins in 
facilitating logic gates between electron qubits^^, and in 
providing ancillas and long-lived memories in large fault- 
tolerant quantum computer architectures. 


State 

Bell Test 

Standard QND 

Tomography 

T PPT C 

ZQC 

2.50(10) 2.70(9) 

97(2)% -0.45(4) 0.88(17) 

DQC 

2.37(12) 2.49(11) 

96(3)% -0.43(6) 0.74(17) 


Table 1 | Summary of electron-nuclear entanglement 
benchmarks. 


METHODS 

The fabrication process, cryogenic setup, gate biasing 
and data acquisition is identical to that described by 
Muhonen et al.^. We combine the signals of three differ¬ 
ent generators (Agilent E8257D/E8267D PSG for ESR 
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and Agilent N5182 MXG for NMR) to apply pulses to 
each of the transitions described in the main text. We 
can achieve ttmwi, rotations in ^ 3 /is and ttrf rotations 
in ^ 30 yus. The generators are pulse modulated using 
a SpinCore PulseBlaster ESR TTL pulse generator. For 
phase control in the density matrix tomography experi¬ 
ments, we use the internal baseband arbitrary waveform 
generator (AWG) in the E8267D vector source, and we 
use an Agilent 81180A AWG to gate the I/Q inputs of 
the N5182. 
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Extended Data Figure 1 | Selective initialization protocol. Top: pulsing scheme used to initialize the 
two-spin system in the state, as described in the main text and Figure Ib-c. The SET current (blue traces) is 
tracked throughout the initialization sequence; high current implies ionized donor. The measurement is discarded if 
the donor is found to be ionized at the end of an electron |^) initialization phase. 
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Extended Data Figure 2 | Bell’s inequality violation experiment details, a. Diagram showing how a ttrf 
maps the parity observable onto the nuclear observable. The pulse coherently swaps the coefficients of the 
|iJ|) ^ \ii[) eigenstates, leaving the odd and even parity manifolds mapped on the |J|) and manifolds 
respectively, b. Correlation measurements used to construct the Bell signal histograms of Figure 2b. Line colours 
correspond to each of the nuclear and electron qubit projection combinations, as shown in the legend {a and (3 are 
defined in Figure 2a). The top inset shows the cumulative moving average of P]y from a sample set of 300 x 4 
single-shot measurements used to obtain one correlation set. For the QND measurements, we calculate the 
correlation E = 1 — 2P^ after each set of 300 single-shot measurements. 
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Extended Data Figure 3 | Density matrix tomography pulse sequences and fitting parameters. Left 
column: set of coherences measured by the tomography protocol. After preparing each tomography target state, a 
ttrf and tt pulse are applied with phase offsets 0 and 0 respectively. These pulses act as phase gates, which alter 
the phase of the coherence. The rest of the pulse sequence maps the coherence to be measured onto the ||J|) 
manifold (column 2). A 7r/2RF pulse is then applied to project the coherence onto the nuclear observable before 
performing a QND nuclear measurement. The nuclear spin proportion is plotted as we increment the phase offsets 
by A0 and A^, and the resulting signals (see Extended Data Figure 4) are fitted to 

= A fp^ph B) C 0.5. Here, Uph is the phase increment number and fp is a frequency in 
cycles/increment which is unique to each coherence (column 3). A and B are free fitting parameters from which we 
extract the off-diagonal element amplitude and phase respectively. The remaining free fitting parameter C is the 
offset of the measured signal, which gives information on the eigenstate populations (column 4). The diagonal 
elements of the density matrix are extracted by constructing a system of equations with all of the extracted offsets 
and I tlT P + I P + I P + I P = 1- The overdetermined system is solved using a non-negative least-squares 
solving algorithm. 
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Extended Data Figure 4 | Density matrix tomography detailed results. Plots correspond to the resulting 
signals obtained from following the density matrix tomography protocol for each of the coherences, after preparing 
the entangled states DQC and ZQC. The data (dots) is fitted (solid line) as described in Extended Data Figure 3. 
The quadrature component is obtained by applying a 7r/2 phase offset to the final projective 7r/2RF pulse and is 
fitted by adding 7r/2 to the sin argument in the fitting function. Note that some traces display oscillations at 
frequencies which are different from the characteristic frequency of the coherence. These oscillations arise from 
cross-talk between coherences that can result from the tomography pulse sequence. The cross-talk is minimized by 
carefully choosing the rate of phase increments for different coherences, to ensure that their Fourier transforms are 
well spaced from each other. Black dashed lines indicate offsets for the ideal states. The resulting density matrices 
for each entangled state are presented in numerical form at the bottom. We assume the matrix is Hermitian and 
apply a global phase correction so the coherence with greatest amplitude has zero phase. 







































